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Abstract. Theoretical spectral energy distributions (SEDs) of white dwarfs provide 
a powerful tool for cross-calibration and sensitivity control of instruments from the far 
infrared to the X-ray energy range. Such SEDs can be calculated from fully metal-line 
blanketed NLTE model-atmospheres that are e.g. computed by the Tubingen NLTE 
Model- Atmosphere Package (TMAP) that has arrived at a high level of sophistication. 
TMAP was successfully employed for the reliable spectral analysis of many hot, com- 
pact post-AGB stars. 

High-quality stellar spectra obtained over a wide energy range establish a data 
base with a large number of spectral lines of many successive ions of different species. 
Their analysis allows to determine effective temperatures, surface gravities, and element 
abundances of individual (pre-)white dwarfs with very small error ranges. We present 
applications of TMAP SEDs for spectral analyses of hot, compact stars in the parameter 
range from (pre-) white dwarfs to neutron stars and demonstrate the improvement of 
flux calibration using white-dwarf SEDs that are e.g. available via registered services 
in the Virtual Observatory. 



1. Introduction 

In an ideal case, the problem of flux calibration can be simplified to 

where fx is the measured stellar flux, A is the radiation area, d is the distance, and Fx 
is the emitted flux. In an astrophysical case, additional absorption in the circumstellar, 
interstellar, and Earth's atmosphere environment will reduce fx- This absorption oc- 
curs on both, time- and energy-dependent scales, and thus exacerbates the calibration 
problem (cf. these proceedings). 

While observational techniques have strongly improved in the last decades, and the 
measurement of fx is performed at high precision, high resolution, and high S/N based 
on ground-based and space-borne instruments, the main uncertainties of Eq.Q]are found 
on its right side. The first challenge to minimize these is to find "simple" objects with 
high photometric stability whose radiation can be described by fundamental physics. 
White dwarfs (WDs) fulfill such requirements. Since their radius R is determined be 
electron degeneracy, their radiative area can be calculated from 
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A = nR 2 (2) 

with 

Ig-m 

(M is the star's mass and g its surface gravity. G is the gravitational constant.)- 

Most of the hot, hydrogen-rich WDs (spectral type DA) with r eff < 40 000 K have 
virtually pure hydrogen atmospheres due to gravitational settling (log g > 7), while the 
hotter WDs exhibit lines of heavier elements due to radiative levitation or weak stellar 
winds. This singled out in the past stars with many metal lines in their spectra as 
standards for flux calibration. 

Moreover, WDs are nearby and, thus, their distances d can be measured with high 
accuracy by parallaxes, particularly with the upcoming mission. A further ad- 

vantage of their small distances is that WDs with extremely low interstellar reddening 
can be found, where the calibration problem reduces almost to the ideal case of Eq.Q]- 
at least for space-borne missions. The spectra of such objects can be used for panchro- 
matic calibration (from the soft X-ray to the far infrared). 

In this paper, we start with a brief summary on the application of theoretical WD 
spectra to the flux calibration (Sect.©. Then (Sect.|3]), we describe modeling of WD at- 
mospheres and the calculation of synthetic spectra. In the case of G 191-B2B is demon- 
strated that state-of-the-art metal-line blanketed model atmospheres reproduce well its 
spectrum (Sect.|4]). In Sect.O we show example for the application of DA model spec- 
tra in the cross-calibration between different instruments, An easy way to access such 
spectra is shown in Sect.[6] We conclude in Sect.[7] 



2. Flux calibration using DA-type white dwarfs 

Most of the presently known (» 10000) WDs are of spectral type DA. Only a few of 
them are relatively bright and exhibit negligible reddening. These are a priori suited to 
act as primary standard stars for the flux calibration in a large variety of observatories. 
As a pre-requisite, however, reliable model atmospheres are necessary to compare ob- 
served with theoretical spectra. These became availab le i n the early 197 Pies, follo wing 
the pioneering work of lAuer & Mihalasl dl969l.ll970l) and lMihalas & Auerl dl970h who 



presented the first NLTE models for p ure hy drogen and hydrogen + helium stellar at- 



ure nydroge i 

mospheres, respectively. Consequently. lOke ( 1974) presented a compilation of absolute 



spectral energy distributions (SEDs) for 38 WDs in the wavelength range A > 3200 A, 
partly up to more than 10000 A. 

Figured] illustrates the advantage of an isolated DA WD in comparison to a DAO 
WD in a binary system - the many strong metal lines in the DAO spectrum hamper 
an accurate calibration in the FUV and the contribution of a cool companion spoils the 
low-energy region. Although the composite synthetic spectrum matches the FUV flux 
as well as the optical and IR magnitudes, the uncertainties in modeling both compo- 
nents cannot be quantified properly. This excludes such objects from the standards of 



' http : //www . rssd . esa . int/index . php?pro j ect=GAIA&page=index 
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Figure 1. Synthetic fluxes of an example binary (A 35, top) and an isolated WD 
(G 191-B2B, bottom) compared with observations. In case of A 35, the composite 
model spectrum is normalized to the 2MASS brightness at 22000 A. The synthetic 
spectrum for G 191-B2B is normalized to that of the calibrated HST spectrum (R. 
Bohlin, priv. comm.). Since the observed and theoretical flux levels are almost 
identical, the observation is plotted in white in the foreground of the TMAP spectrum 
(red) for clarity (visible in the electronic version). The insert shows a sc aled-up part 
of the spectrum. Interstellar reddening was considered using the law of lFitzpatrickl 
(fl999l) . 



panchromatic flux calibration. An isolated DA WD allows to normalize observation 
and synthetic spectrum to the WD brightness in the infrared where interstellar redden- 
ing is lowest. Note that there are no strong interstellar lines in the FUSE spectrum of 
G 191-B2B (Fig.Q}. 
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Oke'slist included HZ43A (WD 1314+293, for the WD numbers cf. lMcCook & Sion 
19991 . and their online catalogue^) and G 191-B2B (WD 0501+527). From a black body 



comparison, Oke found that the observed SED of G 191-B2B corresponds to a black 
body at T = 100 000 K. In addition, HZ 43 A and G 191-B2B were amongst the bright- 
est of Oke's list. T ogether with GD 71 (WD 0549+ 158) and GD 153 (WD 1254+223), 
these were used by lBohlin et all (1 1995b as primary standard stars for the Hubble Space 
Telescope (HST, Tab.[fl. 



Table 1. Parameters of the HST DA standard stars used by iBohlin elaTI d 19951) . 



The otb and ni\ values are from SIMBAD. 



WD name 


name 


?eff 

[K] 


logg 
[cm/sec 2 )] 


ra B 


m\ 


comment 


WD 0549+ 158 


GD71 


32 300 


7.73 


12.78 


13.03 




WD 1254+223 


GD153 


38 500 


7.67 


13.17 


13.40 




WD 0501+527 


G191-B2B 


61300 


7.50 


11.44 


11.69 


metal lines 


WD 1314+293 


HZ43A 


50000 


8.00 


12.52 


12.66 


binary 



The HST flux scale (Bohl in et al.lll995l. 1200 A - 8 000 A) was the reference for 



other calibration work, e.g. Holberg & Bergeroril ( 20061) used it to calibrate four ma- 
jor ground-based systems (Johnson-Kron-Cousins UBVRI, Stromgren uvby, 2MASS0 
JHK S , SDSSQugriz). 

Although stellar model atmospheres improved continuously in the last four deca- 
des, some problems are still unsolved, e.g. the discrepancy between T e s determined 
from Hi lines in ultraviolet and optical spectra dBarstow et al.l 12001 ). In the case of 
G191-B2B. lLaioie & Bergeron! d2007h found = 60680 ± 15 000 K and = 
57414 ± 4 7 00 K. Even improved Stark line-broadening tables for Hi Lyman and 
B aimer lines dTremblay & Bergeronll2009l) leave this discrepancy unexplained. 



The recent, comprehensive analysis of more than 1100 WDs (IGianninas et al. 



201 lh shows significant deviations for the HST DA standard stars (Tab.|2]) from pre- 



vious values (Tab.Q]). Note that the errors given are from their statistical a pproach and, 
thus, appear too optimistic. A calibration of IRAC aboard Spitzer ST dBohlin et al. 



20111) uses a ?W = 61 196 K model for G191-B2B which is within the error limits of 



Gianninas et al.l (1201 lh and fulfils the desired 1 % flux-accuracy limit. 



3. Model atmospheres 

WD spectral analysis requires adequate observatio ns (WDs a re intrinsically faint) and 
state-of-the-art theoretical atmosphere models (cf. IOkdll974h that account for reliable 
physics and deviations from the assumption of a local thermodynamic equilibrium 
(LTE). Figure|2] shows that for cool stars with high s urface gravity (spectral type B 
and later), LTE model atmospheres might be sufficient dAuer & Mihalas|[l972l) . There 



- http : //www . astronomy . villanova . edu/WDCatalog/index . html 
; http : //www . ipac . caltech.edu/2mass/ 
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Table 2. Parameters of the HST DA standard stars (jGiann inas et al. 201 ll) . 



name 


T e ff 

[K] 


logg 
[cm/sec 2 ] 


GD71 


33 590 ± 483 


7.93 + 0.05 


GD153 


40 320 ± 626 


7.93 + 0.05 


G191-B2B 


60920+ 993 


7.55 + 0.05 


HZ43A 


51 116 + 1249 


7.90 ± 0.07 



are, however, NLTE effects in any star, at least towards higher energies and in high- 
resolution spectra. Aiming at a panchromatic flux calibration (at an 1 % flux-error 
limit, which was the designated aim mentioned in many talks at this conference) and a 
cross-correlation of various instruments from the soft X-ray to the infrared, the use of 
NLTE stellar-atmosphere models is mandatory. 




log (7\ ff / K) 



Figure 2. Reliability domains of LTE and NLTE models for static and expanding 
stellar atmospheres. 



Since reliable atomic data is a crucial input for model-atmosphere calculations, 
we start with a brief summary of existing lacks and problems in Sect. 13.11 Then, we 
demonstrate the successful application of NLTE atmosphere models in detailed spectral 
analysis and show recent discoveries (Sect. l3~l2"l) . 
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3.1. Atomic Data 



For TMAP, we compile atomic data fro m several standard s ources like NISTfl Kelly's 



datab ase^), or CHIANTfl the Opacity dSeaton et all 1 19941) and I ron (iHummer et al. 
Il993l) project^] as well as Kurucz's line list^ ( Kurucz 199 ll. 12009 ). The model atoms 
dRauch & Deetienll2003h are provided by the Tubingen Model-Atom Database TMAfJi^ 
(Sect. [6]). TMAD is continuously updated and extended for the most recent atomic data. 
Presently it provides ready-to-use, TMAP-compliant model atoms of the elements H - 
Ca, Ge (Fig.|3]), Kr, and Xe. 




; http : //www .nist . gov/pml/data/asd . cfm 

f http : //www . cf a . harvard . edu/ ampcgi/kelly . pi 

' http : //www . chianti .rl.ac.uk 

f http : //cdsweb . u- strasbg . fr/topbase/home . html 

- http : //kurucz . harvard . edu/ 

" jhttp : //astro . uni - tuebingen . de/ - TMAD 
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3.2. On the reliability of NLTE model atmospheres 



TMAP^ dWerner et alJ2003h can calculate fully metal-line blanketed WD model atmo- 
spheres. About 1600 atomic levels and 7000 individual spec tral lines of the element s 
H - K can be treated in NLTE. Due to a statistical treatment dRauch & Deetjenll2003l) . 
« 300 million of lines of the iron-group elements (Ca - Ni) can be considered. TMAP as- 
sumes hydrostatic and radiative equilibrium and a plane-parallel geometry. It makes use 
of TMAD. TMAP was successfully employed for detailed spectral analyses of hot, com- 



oi i mad. i ma p was, successfully employed to r detailed spectral analyses ot not, com- 
pact stars (e.g. iRauch et all 120071. 12010L (20121 : IWassermann et al.l l2010l : IWerner et al~ 



20121) . In these analyses, the models considered opacities of many elements from H to 
Ni. 

In case of synthetic spectra for high-g DA-type WDs (Sect.Q]), where the calcu- 
lation of pure hydrogen models is sufficient, the highest reliability is achieved because 
H i atomic dat a is accurately known. Especially, tabulated Stark line-broadening tables 
are available (lLemkelll997 : lTremblay & Bergeron! 120091) . The models depend on T e g 
and log g only, i.e. there are no other free parameters. A synthetic WD spectrum can be 
calculated for a large wavelength range with arbitrary resolution. Figure|4]shows the ex- 
ample of EG 274 (WD 1620-391) for the optical to the near infrared wavelength range. 
The observations were performed at ESO (European Southern Observatory, Chile) with 
UVES (Ultraviolet and Visual Echelle Spectrograph) and SINFONI (Sp ectrograph for 
INtegral Field Observations in the Near Infrared) (cf. lVernet et al.|[2008 alfbb . The initial 
aim to achieve an agreement of observation and theory better than 10 % in flux was 
exceeded by far. 




Figure 4. Comparison of a TMAP model (light blue / grey) for EG 274 (r e ff = 
24 276 K, logg = 8.01) with an optical spectrum and infrared brightness measure- 
ments. The insert shows a detail of the UVES spectrum in the optical. The dashed 
lines indicate the flux limits of the synthetic spectrum within the T e g and log g un- 
certainties. 



" |http : //astro . uni- tuebinge n . de/~TMAP| 
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It is worthwhile to note that large observatories like ESO have large calibration 
programs where spectro-photometric standard stars are observed every night. The huge 
number of observations with different setups allows to define master spectra (Fig. @]) 
that help to improve the pipeline reduction of the observations. These spectra may serve 
also to automatically measure e.g. the nights' conditions or to monitor the efficiency 
of a telescope as a whole by comparison of the master flux level of a star to that of an 
actual observation. Moreover, the pipeline reduction can be improved, e.g. to remove 
residuals of echelle orders (see insert in Fig.HJ EG 274 is a spectrum from the Advanced 
Data ProductQ. 

The capability and reliability of a model-atmospheres code can be judged from 
its application to spectral analysis over a wide wavelength range. In general, synthetic 
spectra are more reliable in the low energy range because towards higher energies, the 
consideration of metal opacities (bound-bound as well as bound-free transitions) and 
the lack of atomic input data (Sect. 13.11 ) thwarts any atmosphere code, ho w sophisti 



cated it may be. A simple look at the model-atom stat istics is not sufficient. iKudritzki 



19761) used the complete linearization technique of lAuer & Mihalasl ( 1969 ) and al- 



lowed 18 atomic levels (5 Hi, 1 Hn, 5 Hei, 6 Hen, 1 Hem) to deviate from LTE and 
included six H i line transitions in his so-called classical models. Together with 65 fre- 
quency points, he arrived at a system of 86 equations to solve simultaneously and at the 
accuracy limit of the 32 bit computers that wer e available. Present models employ the 



approximated lambda operator (ALO) method ( Werner & Husfelcfl 1985k IWerne 



1989 1: IWerner & Dreizlerlll999l) and consider all elements from H to Ni ( Rauc 
20031) with e.g. 686 levels in NLTE, 2417 individual lines, and millions of lines of the 



e, 
:h 



1986, 



19971. 



iron group in the analysis of LS V +46°21 dRauch et al.ll2007l) . 

The large number of levels and lines that can be considered in the models, how- 
ever, implies the difficulty to quantify e.g. uncertainties of the respective atomic data 
(Sect. 13.11 ) and, thus, the synthetic spectra. Detailed spectral analyses of high-quality 
spectra and the agreement between model and observation are an important measure to 
estimate the models' flux accuracy. We will now present some individual examples for 
spectral analyses with TMAP models of compact stars with different T e g . 



3.2.1. (Pre-) white dwarfs: T eff » 100 000 K 

The main field of our working group are post-AGB stars, with a focus on hydrogen- 
deficient stellar evolution. The development of observation techniques and the qual- 
ity of the obtained spectra in the ultraviolet wavelength range with e.g. FUSE (Far 
Ultraviolet Spectroscopic Explorer) or HST/STIS (Hubble Space Telescope / Space 
Telescope Imaging Spectrograph) allowed the determination of stellar photospheric pa- 
rameters like ?W, logg , and elemental abundances with unprecedented accuracy (see 
IWerner & Herwigll2006l . and references therein). The identification of hitherto uniden- 
tified lines in FUSE and STIS spectra and abun dance analyses of m any of t he respective 



speci es (e.g. Ne vn (IWe rner et al 



20041) . F vi (IWerner et all 120051) . Arvi dRauch et al 



2007 1 ) . Arvn (IWerner et all 12007 ah . Nevm (IWerner et all l2007bl) Cax dWerner et al 



20081) . and Fe x (IWerner et al.ll2010l). Ga v. G e v, As v, Se iv - vi, Kr vi, vn, Mb vi, Sn iv 



- v, Tevi, I vi, Xevi, vn (IWerner et al.ll2012l) ) were the direct consequence of a con- 
tinuous development of TMAP and the available atomic data. E.g. the identification of 



12 http : //archive . eso . org/eso/eso_archive_adp . html 
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Ne vm lines dWerner et alj|2007b) followed the publica tion of energy levels and spectral 
lines of Ne vm by [Kramida & Buchet-Poulizacl (120061) . 

Figure[5]shows newl y identified Kr and Xe lines in the FUSE spectrum of the DO- 
type WD RE 0503-289 (IWerner et al.ll2012h . With these discoveries, new ionization 
equilibria, e.g. Kr vi / vn, can be evaluated for a more precise determination of T e g. In 
addition, abundance determinations of trans-iron elements establish new challenges for 
stellar evolutionary theory. Although IWerner et al.l ( 2012b identified 22 out of the 23 
Ge v lines that are listed at NIST and five that are listed in Kelly's database, a Ge abun- 
dance anal ysis was impossible due to the lack of transition probabilities. After their 
calculation, iRauch et"al] ( 2012h presented a precise determination of the Ge abundance 
in RE 0503-289 and discovered additional four Ge iv, 14 Ge v and seven Ge vi lines. 



KrVI 
919.934, 



KrVI 
927.34 



Kr VI 
931.39 





Figure 5. Kr vi, Kr vn, Xe vi, and Xe vn lines in the FUV spectrum of RE 0503-289. 
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Figure 6. Two sections of the FUSE (top) and STIS (bo t tom) observations of 
LS V+46°21 compared with the final model of iRauch etal] ( I2007I) . Photospheric 
lines are marked at top of the panels, is denotes interstellar origin (not modeled 
here). 

In the spectral analysis of LS V+46°21 dRauch et alj|2007l) . the central star of 
the planetary nebula Sh2-216, it was possible to identify and reproduce about 95% 
of all spectral lines (more than 1000) in the FUSE and STIS spectra (Fig. [6]) by the 
combination of a photospheric model and a model of the interstellar line-absorption 
Owens profile fitting procedure. 

Figure [6] demonstrates the unique opportunity to employ stars as stellar laborato- 
ries for atomic physics. The large number of spectral lines, e.g. of Fe vi, allows to 
measure relative line strengths and to verify or improve the respective log gf values 
(g is here the statistical weight of the lower level of the line transition, / the oscilla- 
tor strength). In addition, precise wavelength measurements in high-resolution spectra 
allow to improve the atomic level energies. 



3.2.2. Super-soft sources: T e « x 700 000 K 

Some months after their outbursts, novae may evolve into very bright super-soft X-ray 
sourc es (SSS). V4743 Sgr (outburst Sep 20, 2002) was the brightest SSS in March 



2003 (INess et al.ll2003l). Although th e origin and the physics of the nova's spectrum is 



still under debate, Rauch et all (1201 Ol) could reproduce well the strengths and shapes of 



lines and absorption edges and determined r e ft- = 740000 ± 20 000 K. Fig. [7] shows a 
comparison of theoretical spectra with an XMM-Newton (X-ray Multi-Mirror Mission) 
RGS (Reflection Grating Spectrometer) fluxed spectrum from BiRE^. Note that the 
given error is rather conservative, i.e. we arrive at an uncertainty level of better than 
3 % even for these extremely hot objects. 



13 http : / /xitim . esac . esa . int/BiRD/ 




Figure 7. Dependance of the strengths of lines (marked above the top spectrum, 
red) and ground-state thresholds (marked at top) on T e ff . The SEDs of our models 
(log g = 9) are compared with the XMM-Newton observation (2003-04-04) of V 4743 
Sgr and normalized to match the observed flux level at 38 A. IG denotes a generic 
model atom that contains Ca - Ni. 



V 4743 Sgr is not the only example where our static models were used success- 
fully to determine photospheric parameters of the WD in SSS. Nova KT Eri (outburst at 
Nov 25, 2009) was observed in early 2010 with Chandra and XMM-Newton. In Figj8j 
we show a comparison of a model SED with the observation. Although the model is 
not fine-tuned in abundances and does not include all elements from H to Ni, the re- 
production of the observed spectrum is already very good. Obviously, lines of different 
ionization stages, e.g. N vi and N vu are formed in different radial velocity regimes but 
their theoretical strengths match the observation. 



3.2.3. Neutron stars: T eff s 10 000 000 K 

At a typical log g * 14.39 (neutron star with M = 1.4M and R = 10 km), an extended 
Fig. H would show that TMAP models are still valid up to 7/ eff ~ 20 MK. At higher r eff , 
comptonization has a significant impact on the X-ray spectrum ( Rauch et al.1 20081). 

In the case of EXO 0748-676, an X-ray bursting neutron star, iRauch et al.l (120081) 
could easily show that the identifica t ion of gravitationally red-shifted absorption lines 
of Fe xxv and xxvi by ICottam et all (120021) is wrong because at the relevant tempera- 
tures, Fexxiv lines should be much stronger but these were not observed (Fig. [9]). 
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Figure 8. SED of our preliminary TMAP model (logg = 9) compared with the 
Chandra spectrum (2010-01 -31) of KTEri (normalized to the observed flux at 38 A). 




10 15 20 25 

wavelength / A 



Figure 9. 



Model spectra for EXO 0748-676 with different iron content. 
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4. G 191-B2B - too many metal lines? 



G 191-B2B is the visually brightest of the HST standard stars (Table© and thus it may 
be the most observed and best studied hot DA WD in all wavelength ranges. Hence, it 
should be the prime standar d star for the UV wa velengths as well. However, it shows 
many weak metal lines (e.g. iBarstow et al.ll2003T) in its observed UV spectrum. There- 
fore its was not used for the FUSE calibration (G. Sonneborn, this conference). De- 
tailed modeling with TMAP shows that these lines can be well reproduced and, thus, a 
reliable synthetic spectrum from the H i Lyman threshold to the far IR can be calculated 
that allows a flux calibration over a wide energy interval. Fig. Q] shows a comparison of 
HST observations of G 191-B2B0 within 1 400 < A < 30000 A with our preliminary 
TMAP model that agree almost perfectly in the continuum-flux level. A detailed fine- 
tuning of the metal abundances is still ongoing that will improve the agreement in the 
FUV and EUV range. 

Figure [TU] shows a comparison of the flux of a pure hydrogen model to a metal- 
blanketed model and a black-body flux. Strong metal-line blanketing is visible around 
the flux maximum in the H-Ni model. This results in significant deviations from the 
pure H model flux even in the IR. Compared to the model-atmosphere fluxes, the black 
body has its flux maximum at a 2.7 times higher wavelength and the peak maximum is 
about 40% lower. 




Figure 10. Astrophysical flux of two TMAP models for G191-B2B with 
different elemental compositions (blue, thick: pure H atmosphere; red, thin: 
H+He+C+N+O+Si+P+S+Ca+Sc+Ti+V+Cr+Mn+Fe+Co+Ni). The dashed line 
shows the black body flux at the same T = 60 920 K. 



H http : //www. stsci . edu/hst/observatory/cdbs/calspec . html 
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5. Cross-correlation between various instruments 



Beuermann et all (120061) established two DA WDs, namely HZ43 A and SiriusB, and 
the neutron star RXJ1 85635-3754 as soft X-ray standards. This enabled a cross- 
calibration between the Chandra LETG+HRC-S (Low Energy Transmission Grating, 
High Resolution Camera Spectroscopic array), the EUVE (Extreme Ultraviolet Ex- 
plorer) spectrometer, and the ROSAT PSPC (Rontgensatellit, Position Sensitive Pro- 
portional Counter). With their attempt, the effective area of ROSAT could be corrected, 
yielding consistent fluxes from about 20 to 160 A. It is worthwhile to note that dur- 
ing the work for this cross-correlation, Jelle Kaastra discovered discrepancies between 
TMAP and TLUSTY dHubeny & Land 1 19951) mo dels in the ca se of HZ43 A, that were 
then intensively in vestigated and understood (e.g. lRauchll2008l) . 

brake! (120121) and the Chandra LETG team relied on the growing confidence in 
the photospheric parameters for HZ43 A dBeuermann et al.|[2006l) and that its spectrum 
was well-represented by pure hydrogen atmospheric models. They used HZ43 A spec- 
tra obtained over ten years and TMAP atmospheres in order to improve the in-flight 
calibration. The corrected effective area curve allowed then to incorporate a secular 
decline in quantum efficiency for a better (re-)calibr ation of LET G spectra. It is worth- 
while to note that GD 153 was newly introduced by lMenzl ( 2011 ) for the calibration of 
LETG spectra in order to bridge the T e g gap between Sirius B and HZ43 A (Tabled. 



6. Virtual Observatory Services and Tools for Synthetic Spectra 

In the framework of the Virtual ObservatorJPl (VO), the German Astrophysical Virtual 
Observator¥^ (GAVO) provides the registered VO service TheoSSA (Theoretical Spec- 
tra Simple Access, Fig.fTTb. It was created to provide access to synthetic stellar SEDs at 
three levels. 1) Fast and easy download of precalculated SEDs by an unexperienced VO 
user (no detailed knowledge of the model-atmosphere code necessary). 2) Calculation 
of SEDs (elements H, He, C, N, and O) with standard TMAD model atoms for a prelim- 
inary analysis of individual objects via the TMAP web interface TMAWO 3) Advanced 
calculation of SEDs with especially tailored model atoms, e.g. for a comparison with 
other stellar model-atmosphere codes. 

TheoSSA is related to the GAVO database that includes already thousands of SEDs. 
Newly calculated SEDs are automatically ingested and, thus, this database is growing 
in time. If a requested SED is not available, the VO user will be redirected to the TMAW 
interface, where individual models can be calculated. Figure[[2]illustrates the dataflow 
between the TheoSSA and TMAW WWW interfaces and the GAVO database. For a pure 
H model, the calculation takes some hours, for H+He+C+N+O composed models, one 
or two days are necessary. The wall-clock time, however, is depending on the number 
of queued SED requests. The calculation of extended model grids (> 150 models) 
makes use of compute resources of AstroGrid-dB- 
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" |http : //www, g-vo" . org/ 

1 ]http : //astro . uni - tuebingen . de/ - TMAW 

" jhttp : //www. gac- grid . de/ 
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G AVO TheoSSA Web Interface 



Help 



Service info 
Related 

Compute custom 
SEDs 

Metadata 

ldertrfier» 
Description » 
Keywords » 
Creator » 
Created » 
Data up dated » 
Copyright» 
Source » 
Reference URL» 



Ti-y AOX to query 
□ ur data. 



TheoSSA provides spectral energy distributions based on model 
atmosphere calculations. Currently, we serve results obtained using 
the Tubingen NLTE Model Atmosphere Package for hot compact 
stars. 



- = = :: report errors 



Effective between and 

Temperature Rartge of tiie atmosphefe-'s effects temp era Cures to include. If you only specify one bound, you get a half-infinite internal. 
[K\ 

Log. Surface between and 

gravity [cm/s2] Range of surface gravities to include. If you only specify one bound, you get a half-infinite interval. 



Mass loss rate between and 

[Msun/yr] Range of mass loss rates to include. If you only specify one bound, you get a half-infinite interval 



Mass Fraction "any t " between and 

1 Mass fraction of an element. You may leave out either upper or lower bound. 



Mass Fraction ant t between and 

2 Mass fraction of an element. You may leave out either upper or tower sound. 



Mass Fraction ant t between and 

3 Mass fraction of an element. You may 'esve oa! either upper or lower bound. 



Standard Stars gisi-es - 

Common name of object observed. 

Type » VOIaWe 

O Text 

Format of spectrum data files 



Table Sort by T Limrtto im T items. 

Output format html | Mramgwirafc | 



Result Imkl £ 



Figure 11. TheoSSA web interface at http : //dc . g-vo . org/theossa 



C 



^TheoSSA request: T e ff, log g, 

I - 

check GAVO database: 

requested parameters matched 
within tolerance limits? 
I 




yes: offer existing SED 



X 



no: calculate new model and SED 
and ingest into GAVO database 



request exact parameters 




GAVO database 



^eif, log g, X, 
metadata 



GAVO database 

r cff , logg,Xj 

models 
atomic data 
frequency grids 
SEDs 



Figure 12. TheoSSA / TMAW data flow. X\ is the mass fraction of species i e 
[H,He,C,N,0]. 



Presently, TheoSSA provides a preliminary sample of pre-calculated DA WD spec- 
tra (Table [3]) from the optical to infrared wavelength range. These will be expanded to 
include the X-ray and UV ranges in the near future. 
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Table 3. 


Standard star SEDs (3 000 


- 55 000 A) presently available in TheoSSA. 


name 


1 aff 

* en 


l02 8 

rprn / QPr*^l 


comments 




24276 


8.01 




Feige 67 


1C AAA 

/j uuu 






ppi pp 1 i n 

1 ^±£i^_. 1 1U 


40000 


5.00 


He mass fraction 107 


GD108 


22 908 


5.30 




GD153 


38 205 


7.89 




GD50 


40 550 


9.22 




GD71 


32 780 


7.83 




G191-B2B 


60920 


7.55 


metal-line blanketed model in preparation 


G 93-48 


18100 


7.85 




HZ 2 


20 600 


7.90 




HZ 43 


51 116 


7.90 




Sirius B 


24 826 


8.60 





7. Conclusions 

DA-type WDs are stable sources of radiation on a wide energy range from the X-ray to 
the IR. This allows a cross-correlation between ground- and space-based observations. 

WD modeling arrived at a high level of sophistication and thus, they are well 
suited to provide model SEDs within 1 % flux-level accuracy. It is, however, difficult to 
quantify the absolute accuracy of model atmospheres, at least if we aim at the ambitious 
1 %. About 200000000 variables are calculated iteratively to achieve convergence of 
a NLTE stellar-atmosphere model. Besides the mentioned uncertainties in atomic data 
and line-broadening theory, systematic deviations may exist between different codes 
due to programming language, coding, compiling, properties of the compute nodes, 
etc. Future close collaboration of all group s working in stellar-atmosphere modeling 
is required (in the legacy of past workshops. ICrivellari et al.|[l99ll : lHubeny et alj|2003 ) 



in order to identify and minimize existing uncertainties and to further improve the flux 
calibration which is a continuous process. The VO provides an ideal platform to publish 
observations as well as models. The quality of spectral analysis and, hence, stellar- 
atmosphere models benefits from the observational data that it is based on. 

Further improvement of flux calibration will be possible with GAIA. Precise dis- 
tances of almost all known WDs will be available in the post-GAIA era. 

For ground-based telescopes, it is a major challenge to model and to correct for 
Earth's atmosphere contributions (cf. a variety of papers in these proceedings) and to 
establish an extended database of secondary standard stars, that should be observable 
in the vicinity of science targets. 
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